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ABSTRACT 



Aims. We present a study of the broad optical/UV emission line profiles in AGN (active galactic nuclei) to get information on the 
dominant velocity components (turbulence, rotation, etc.) in the central broad-line region (BLR). 

Methods. We introduce line broadening simulations of emission line profiles and compare these results with the largest homogeneous 
data set of reverberation-mapped AGN. 

Results. The underlying broad-line profiles in AGN are Lorentzian profiles caused by turbulence in the line emitting region. The 
turbulent velocities are different for the different line emitting regions of Hy, Ha, Lya, Cm]/I1909, Heii/11640, and SiiV/11400. 
The turbulent velocities go from 400 km s^ 1 for H/J up to 3,800 km s~' for Lya+N v A 1240. The dominating broadening mechanism 
of these profiles is broadening due to rotation. The rotation velocities causing the line profile broadening go from 500 km s~ l up to 
6,500 km s~' . Here we present interrelations between observed emission line widths (FWHM) and their related rotational velocities to 
correct for the contribution of the turbulence to the broad-line profiles. 
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1. Introduction 

It is now generally accepted that a super-massive black hole in 
the center of active galactic nuclei (AGN) is surrounded by an 
accretion disk (e.g. Urry & Padovani 1995). In the outer region- 
sof the accretion disk, the broad emission lines are created by 
photoionization. However, many details of this line emitting re- 
gion are unknown, and there are many models that treat their ge- 
ometry and structure (e.g. Collin-Souffrin et al. 1988 Emmering 
etal. [19921 Konigl & Kartj e [T994l M urray & ChiangfTWTl Elvis 
120001 Kollatschny I2003al Ho [20081 Gaskell l20T0l Goad et al. 
120 121 and references therein). This broad-line region is spatially 
not resolved on direct images. 

Information on the structure and kinematics of this region 
can in principle be obtained from the broad emission line pro- 
files. However, there is the ambiguity in the possibility that dif- 
ferent geometries and different kinematics result in emission line 
profiles that have similar shapes. This problem makes it difficult 
to disentangle the profile contribution of individual parameters 
(e.g. Netzer [19901 Sulentic et al. 120001 Zamfir et al. [20101 and 
references therein). The shape and the width of the emission line 
profiles in AGN might depend, among others parameters, on the 
velocity field, on the geometry of the line emitting gas, on obscu- 
ration effects, on the superposition of line emission from differ- 
ent regions, and on the isotropy/anisotropy of the line emission. 

The velocity field might be a superposition of different 
components, such as Doppler motions, turbulence, shock com- 
ponents, in/outflow components, and rotation. Different veloc- 
ity components result in different profiles, and the final pro- 
file is a convolution of different components. Doppler mo- 
tions in the gravitational potential of a central black hole can 
result in Gaussian profiles, where the width of the lines re- 
flects the Doppler motion of the line emitting gas (Peterson 
& Wa ndel [19991 Fromerth & Melia 120001 Peterson & Wandel, 
12000b . Logarithmic profiles agree well with observations of 



emission line wings and can be produced by outflowing gas (e.g. 
Blumenthal & Mathews [19751 Netzer [T990l l. The existence of 
Lorentzian profiles is consistent with emission from extended 
accretion disks (Veron-Cetty et al. 120011 Sulentic et al. I20021 >. 
Goad et al. ( 120121 1 demonstrate that turbulent motions in the 
outer accretion disk produce Lorentzian profiles. Furthermore, 
Lorentzian profiles might be explained by shock breakouts in 
dense winds (Moriya & Tominaga 12012b . In addition to these 
base-line profiles, both rotation and electron scattering (e.g. Ha 
in NGC 4593; Laor 2006) can lead to broadening of the line pro- 
files. 

Osterbrock pointed out as early as 1 19781 that a combination 
of turbulence with rotation agrees well with the distribution of 
observed line widths of Seyfert 1 galaxies known at that time. 
More recently, researchers have explored the possibility of an 
additional net infall component (Gaskell 120101 for a review, Hu 
et al. 120 121 but see Sulentic et al. !20121 > and an additional outflow 
component in luminous AGNs (e.g. Wang et al. 1201 11 1. 

The contribution of the velocity components might be differ- 
ent from line to line, as well as from galaxy to galaxy. In their 
study of AGN line profiles Robinson et al. (1995 ) demonstrated 
that different types of models produce similar kinds of profiles 
and that the profile shapes differ between the AGN populations. 
In addition, the optical and UV emission line profiles are simi- 
larly diverse in appearance. In an investigation of broad emission 
lines from the Sloan Digital Sky Survey (SDSS) Richards et al. 
( 2002) point out that the explanation of their profiles is one of the 
most important tasks of future AGN models. Using CIV emis- 
sion line profiles, Richards et al. (1201 11 1 show that the differences 
in the accretion disk wind between quasars can account for some 
of the diversity of broad emission line profiles. 

We demonstrated in a first paper (Kollatschny & Zetzl 120 lTl 
hereafter called Paper I) that there is a general trend between 
the full-width at half maximum (FWHM) and the line width ra- 
tio FWHM/cr;„, e in the broad emission lines of AGN. Different 
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emission lines exhibit different systematics in the FWHM/o"ji„ e - 
vs-FWHM diagram. The line width FWHM reflects the rota- 
tional motion of the broad-line gas in combination with the as- 
sociated turbulent motion. This turbulent velocity is different for 
the different emission lines. 

Here we present a further study of additional broad optical 
and UV emission line profiles of reverberation-mapped AGN. 
The line profiles are parameterized by the ratio of their full- width 
at half maximum (FWHM) to their line dispersion au ne . 

2. Data: The line profile sample 

We studied in detail the mean profiles, as well as the root-mean- 
square (rms) profiles of the H/3, He n ,14686, and C iv A 1550 lines 
of an AGN sample in Paper I. This set of observations is the 
largest homogeneous sample of variable AGN at this time, and it 
is based on reverberation-mapped AGN spectra (Peterson et al., 
I2004t . The optical spectra were collected with different ground 
based telescopes while the UV spectra were taken with the IUE 
satellite and the Hubble Space Telescope. The sample has the 
advantage that all the spectra were reduced in exactly the same 
way and that all spectra of each galaxy were inter-calibrated with 
respect to each other. The original AGN sample consists of 37 
objects. In some cases multiple line spectra exist for a particu- 
lar galaxy based on different variability campaigns of the same 
galaxy. 

The present study is based on this data set as well. Here we 
analyze all the remaining optical and UV broad-line profiles we 
did not consider in Paper I, i.e., the Hy, Ho-, Lya, Heii/11640, 
C m] ,11909, and S iv ,11400 lines. 

Usually the observed broad-line profiles in AGN are more 
or less contaminated by additional narrow emission line compo- 
nents from the narrow line region. For avoiding any major ambi- 
guity we only inspected the root-mean-square (rms) line profiles 
out of the sample of Peterson (120041 1 as we did in Paper I. The 
rms profiles display the clean profiles of the variable broad emis- 
sion lines. The narrow line components disappear in these spec- 
tra because they are constant on time scales of years. The optical 
and UV emission line profiles of our AGN sample are parame- 
terized by their line widths (FWHM, and cr;„, e ). The relationship 
between FWHM and cr/ lne contains information on the shape of 
the profile. We use all those line profiles of the sample that were 
regarded as both reliable and less so by Peterson d2004t - same 
as in Paper I. 

All spectral data we are using in this current investigation 
are listed in Table 1. We present Hy, Ha, Lya, Cm] A 1909, 
Si iv ,U400, and He n ,11640 line widths (FWHM) and line width 
ratios (FWHM/cr) of all galaxies in the AGN sample. 

3. Results 

The following are the emission line profiles resulting from vari- 
ous kinematical and dynamical models that have been discussed 
in the context of the BLR in AGNs: 

- Gaussian profiles due to Doppler motions: 

FWHM/cr,,„ e = 2.35; 

- Lorentzian profiles due to turbulent motions: 

FWHM/cr i ine * 1 (-» 0); 

- Exponential profiles caused by electron scattering: 

FWHM/cr,,„ e = V21n2 * 0.98; 

- Logarithmic profiles caused by in-/outflow motions: 

FWHM/cr fe * 1 (-» 0). 
We present in Fig. 1 all these emission line profiles scaled to 



the same H/3 line width (FWHM). The exact mathematical value 
of FWHM/o7,„ f goes to zero for Lorentzian profiles because of 
the infinitely extended wings (see Fig. 1). However, the typ- 
ical FWHM/cr;„, e values we are measuring are on the order 
of one when integrating the spectral lines over a few hundred 
Angstroms. Similarly the exact FWHM/cr /( „ e value for logarith- 
mic profiles cannot be calculated because the central intensity 
goes to infinity. See also discussions of this ratio in Peterson et 
al. d2004l l. Collin et al. d2006l and Goad et al. d20T2l l. 

In Fig. 2 we again show the basic logarithmic H/3 profile for 
Vminlvmax = 0.003, as well as more line profiles generated by ex- 
panding spherical shells of radiatively-driven isotropically emit- 
ting clouds with different v mi „ /v max values ranging from 0.01 up 
to 0.9999 (see also Capriotti et al . [19501 ■ 

3.1. Line profile broadening simulations 

We are trying to find appropriate answers to the following ques- 
tions with our present investigation of AGN emission line pro- 
files: which line profile is emitted intrinsically and what broad- 
ening mechanisms (besides the instrumental broadening) have 
an impact on the observed profile. Lorentzian emission line pro- 
files and/or Gaussian profiles are the most accepted profiles that 
are thought to be emitted intrinsically. 

We showed in Paper I that rotational line broadening is in- 
deed the most important broadening parameter for AGN emis- 
sion line profiles. We calculated the rotational line broadening 
by the convolution of Lorentzian or Gaussian profiles with ellip- 
soidal profiles. The rotational velocity b = AA/x is by definition 
the half width at zero intensity (HWZI) of an ellipsoidal profile 

A(x) = - Vl - x 2 . (1) 
n 

In Paper I we used a numerical code that was developed 
by Hubeny et al. (119941 1 to compute the line broadening due to 
rotation. This program considers limb darkening as well. In the 
present study we use our own routine without limb darkening 
that is based on the following broadening formula: 

S(y) = \ W(y- x)A(x)dx (2) 



(Unsoeld. 119551 , where W is the intrinsic line profile without ro- 
tational broadening, A the rotational profile, and S the convolved 
profile. 

Other intrinsic line profiles and other broadening mecha- 
nisms than rotation lead to different profile shapes. We computed 
these profile shapes as well and compared them with observed 
AGN profiles. 

Figures 3 and 4 show computed Lorentzian profiles that were 
broadened by rotational motions. The rotation velocities range 
from 100 km s' 1 to 8,000 km s . For the intrinsic H/3 Lorentzian 
profile we adopted a turbulent velocity of 500 km (see Paper 
I) in the line emitting region (Fig. 3). For the intrinsic Lorentzian 
profile of the C iv A 1549 line we adopted a turbulent velocity of 
the line emitting region of 3,000 km s (Fig. 4). We present ad- 
ditional line broadened profiles in Figs. 5 and 6. In Fig. 5 we 
made the assumption that the intrinsic H/3 profile is a Gaussian 
profile (FWHM = 500 km s^ 1 ) that is broadened by rotation 
as in Fig. 3. Figure 6 demonstrates the line broadening of a 
Lorentzian H/3 profile owing to Doppler motions (convolution 
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Table 1. Hy, Ha, Lyo-, C m] ,11909, Si iv ,11400, and He n ,11640 line widths of our AGN sample. 
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Object 


Line 


FWHM 


0" 




FWHM/tr 






[km s~ 


'] 


[km s~ 


'] 






PG 1229+204 


Hy 


3207 ± 


875 


1540 ± 


159 


2.082 


±0.169 


PG0844+349 


Hy 


4946 ± 


1085 


2058 ± 


218 


2.403 


±0.195 


PG0052+251 


Hy 


5633 ± 


3585 


2230 ± 


502 


2.526 


±0.432 


PG0953+414 


Hy 


2960 ± 


589 


1299 ± 


193 


2.279 


±0.261 


PG1307+085 


Hy 


4278 ± 


881 


1758 ± 


193 


2.433 


± 0.204 


PG2 130+099 


Hy 


2661 ± 


481 


1836 ± 


191 


1.449 


±0.130 


PG 1226+023 


Hy 


3274 ± 


484 


1688 ± 


142 


1.940 


±0.130 


PG0804+761: 


Hy 


2201 ± 


295 


1256 ± 


141 


1.752 


±0.160 


PG1211 + 143: 


Hy 


2014 ± 


249 


1376± 


157 


1.464 


±0.143 


PG0844+349 


Ha 


2436 ± 


329 


1625 ± 


73 


1.499 


± 0.057 


PG1426+015 


Ha 


5450 ± 


842 


4254 ± 


290 


1.281 


±0.078 


PG1229+204 


Ha 


3229 ± 


364 


1737 ± 


118 


1.859 


±0.101 


PG1617+175 


Ha 


3794 ± 


780 


2483 ± 


160 


1.528 


±0.083 


PG 14 11+442 


Ha 


1877 ± 


375 


2437 + 


196 


0.770 


± 0.072 


PG0026+129 


Ha 


1117 ± 


109 


1961 ± 


135 


0.570 


±0.056 


PG0804+761 


Ha 


3155 ± 


569 


2046 ± 


138 


1.542 


±0.088 


PG2 130+099 


Ha 


1574 ± 


438 


1421 ± 


80 


1.108 


±0.059 


NGC4151 


Ha 


3156 + 


300 


2422 ± 


79 


1.303 


±0.038 


NGC4593 


Ha 


3399 ± 


196 


1253 + 


90 


2.713 


±0.147 


NGC7469 


Ha 


1615 ± 


119 


1164 + 


68 


1.387 


±0.071 


NGC3516 


Ha 


4770 ± 


893 


2108 + 


69 


2.263 


±0.057 


NGC3227 


Ha 


3168 ± 


67 


1977 + 


134 


1.602 


±0.091 


PG0052+251: 


Ha 


2682 ± 


453 


1913 + 


85 


1.402 


±0.054 


NGC4151: 


Ha 


3724 ± 


529 


1721 + 


47 


2.164 


±0.046 


PG1211+143: 


Ha 


1425 + 


382 


2321 + 


231 


0.614 


±0.083 


PG 1226+023: 


Ha 


1638 ± 


424 


2075 + 


239 


0.789 


±0.104 


PG1307+085: 


Ha 


3084 ± 


1041 


1843 + 


98 


1.673 


±0.073 


Mrk279: 


Ha 


3408 ± 


555 


1405 + 


266 


2.426 


±0.351 


NGC5548: 


Ha 


3044 ± 


381 


1694 + 


80 


1.797 


±0.069 


Fairall9 


Lya 


3503 ± 


1474 


4120 + 


308 


0.850 


±0.069 


3C390.3 


Lya 


8732 ± 


985 


3952 + 


203 


2.210 


±0.088 


3C390.3: 


Lya 


8225 ± 


781 


4600 + 


141 


1.788 


± 0.044 


NGC5548 


Cm] ,11909 


4895 ± 


1263 


3227 + 


176 


1.517 


± 0.070 


NGC5548 


Cm] ,11909 


5018 + 


1458 


2360 + 


222 


2.126 


±0.156 


NGC7469 


Si iv A 1400 


6033 ± 


1112 


3495 + 


269 


1.726 


±0.109 


NGC3783 


Si iv ,11400 


6343 ± 


2021 


3488 + 


161 


1.819 


±0.068 


NGC5548 


Si iv ,11400 


7044 ± 


1849 


4014 + 


253 


1.755 


±0.090 


NGC5548 


Si iv ,11400 


6455 ± 


3030 


2576 + 


389 


2.506 


±0.288 


NGC7469 


He nil 640 


10725 ± 


1697 


3723 + 


113 


2.881 


±0.065 


NGC3783 


He n A 1640 


8008 ± 


1268 


3870 + 


162 


2.069 


±0.068 


NGC5548 


He n A 1640 


8929 ± 


1571 


4397 + 


154 


2.031 


±0.056 


NGC5548 


He ii ,11 640 


9803 ± 


1594 


3897 + 


264 


2.516 


±0.130 



Note. Measurements of objects with colons (:) were considered as less reliable by the authors (Peterson et al.. 12004b . 



with a Gaussian profile). In all cases we made the assumption 
that the integrated line intensities remain constant. 

3.2. Observed and modeled line width ratios FWHM/ au ne 
versus line width FWHM. 

In the next step we present a comparison of the observed line 
profile shapes in different AGN spectra with computed line pro- 
files. 



3.2.1. Rotational line broadening of Lorentzian H/3, 
He ii ,14686, C iv ,11549 profiles 

In Paper I (their Figs. 1-3) we showed the observed H/?, 
Hen ,14686, and Civ Al 549 line width ratios FWHM/cr/,„ e ver- 
sus their line width FWHM in AGN. These line width data were 
taken from the sample of Peterson et al. ( 120041 ). We likewise in- 



cluded those measurements that were regarded as less reliable 
by the authors because they follow exactly the same trend. 

We modeled the observations by rotationally broadened 
Lorentzian profiles. We present in Fig. 7 the observed ratios 
FWHM/cr,,„ e vs. FWHM for the H/3, He n /14686, and C iv ,11549 
profiles in one single plot, as well as their modeling with our 
new routine. The rotational velocities in this figure are slightly 
shifted towards higher FWHM (in comparison to Paper I) be- 
cause the effect of limb darkening leads to a slightly modified 
rotational profile (see Fig. 168 in Unsoeld, [1955| . The data for 
these three line profiles and their corresponding models fill dif- 
ferent areas in this plot. The underlying H/3 Lorentzian profile 
has a line width of 4-00 km s~ l , the Hen ,14686 has an intrinsic 
line width of 900 km s~ x , and the C iv /1 1549 line has an intrinsic 
line width of 2,900 km s . These line widths are least square fits 
to the data in Fig. 7. 

The modeled rotational velocities of the H/3 lines range from 
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Velocity [kms ] 

-1500 -1250 -1000 -750 -500 -250 250 500 750 1000 1250 1500 



Velocity [kms 1 ] 




Gauss (Doppler) 
Lorentz (turbulence) ■ 
Exp (e~ scattering) 
lion (irWoutflow) 
(v mjn /v max =0.003) 



Wavelength X [A] 

Fig. 1. Basic emission line profiles resulting from different 
kinematical models for the BLR in AGN: Gaussian (solid line, 
blue), Lorentzian (dashed, black), exponential (dot dot dashed, 
red), and logarithmic (dot dashed, green) profiles. All profiles 
are scaled to the same H/3 line width (FWHM=500 km s~ l ). 




Hp FWHM v log =500 kms -1 



,=0.003 

=0.01 

=0.05 

-0.1 

-0.2 

.0.3 

.0.4 

.0.5 

=0.6 

-0.7 

=0.8 

=0.9 

.0.9999 



4850 4900 
Wavelength X (A] 



Fig. 2. Logarithmic H/3 profile with an intrinsic line width 
(FWHM) corresponding to 500 kms -1 . Furthermore, we cal- 
culated the emission line profiles caused by expanding spheri- 
cal shells of radiatively-drivenisotropically emitting clouds for 
different values of v„„„ to v max (see also Capriotti et al. |1980l l. 



Velocity [kms 1 ] 
-2500 2500 





Hp FWHM v^^OO kms" 1 




+v rDl = 1 00 kms"' 






= 500 kms" 






=1000 kms 






=2000 kms 






=3000 kms 






=4000 kms _ 






=5000 kms — 






=6000 kms 






=7000 kms 






=8000 kms 

L ■ 





4850 
Wavelength X [A] 



Fig. 3. Line broadening of a Lorentzian H/3 profile due to 
rotation. The rotation velocities range from 100 km s^ 1 to 
8,000 kms~ l . The intrinsic turbulent velocity of the H/3 line 
corresponds to 500 km s -1 . All profiles are normalized to have 
the same total flux. 



5000 7500 10000 
CIV FWHM vl llr+) =3000 kms -1 — 




1540 1560 
Wavelength X [A] 



Fig. 4. Line broadening of a Lorentzian C iv /1 1549 profile due 
to rotation. The rotation velocities range from 100 Aim s^ 1 up to 
8,000 kms~ l . The intrinsic turbulent velocity of the C iv A 1549 
line corresponds to 3,000 kms~ l . 




Velocity [kms 1 ] 
-2500 2500 



4850 4900 
Wavelength X[A] 



Fig. 5. Line broadening of a Gaussian Wfi profile due to ro- 
tation. The rotation velocities range from 100 km s~ l up to 
8,000 fan J -1 . The intrinsic line width of the H/3 line corre- 
sponds to 500 km s~ l . 




4850 4900 
Wavelength X [A] 



Fig. 6. Line broadening of a Lorentzian H/3 profile due to 
Doppler motions. The intrinsic line width of the H/? line cor- 
responds to 500 A77? s . This Lorentzian profile has been con- 
volved with Gaussian profiles having widths of 100 few s~ l up 
to 8,000 km s~ l . 
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5 



3 




12000 



FWHM [kms~ 



Fig. 7. Observed line width ratios FWHM/ cri me vs. line width 
FWHM for H/3 (black diamonds), Hen ,14686 (blue squares), 
and Civ/11549 (red circles). The dashed curves represent 
their corresponding theoretical line width ratios based on ro- 
tational line broadened Lorentzian profiles (FWHM = 400, 
900, and 2,900 km s~ l ). The rotation velocities go from 500 to 
7000 km s~ l (curved dotted lines from left to right). A color ver- 
sion of this figure is available in the online journal. 
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x 
3 



3000 




5000 



10000 



12000 



FWHM [kms~ 



Fig. 8. Observed H/3 line width ratios FWHM/ cr ftne vs. line 
width FWHM. The curves represent their modeling by rota- 
tional broadening (500 - 7,000 kms~ l ) of Gaussian profiles 
(500, 1,000, 2,000, 3,000, 4,000, 5,000km s~ l line width). A 
color version of this figure is available in the online journal. 
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FWHMv )lirh 50 kms~ 
100 
200 
500 
1000 
2000 



2000 



4000 6000 8000 10000 12000 



FWHM [kms~ 



Fig. 9. Observed H/3 line width ratios FWHM/ crime vs. line 
width FWHM. The curves represent their modeling by convolv- 
ing Lorentzian profiles (50, 100, 200, 500, 1000, 2000 km s' 1 
line width) with Gaussian profiles (500 - 1 0,000 kms~ l line 
width). A color version of this figure is available in the online 
journal. 



c 

5 
I 




12000 



FWHM [kms~ 



Fig. 10. Observed H/3 line width ratios FWHM/ cri ine vs. line 
width FWHM. The curves represent their modeling by logarith- 
mic profiles for different v m ,„ / v max ratios (0.01 - 0.9, from bot- 
tom to top). Furthermore, we varied the line width from 300, 
500, 700, 900, to 1 100 km sT 1 (left to right). 



500 to 6,000 km s in the individual galaxies. The 
He n /14686 and C iv /11549 lines exhibit higher rotational veloc- 
ities of at least 1,500 km s~ l . 

3.2.2. Fitting the observed H/3 line width ratios with other 
models 

We generated further test models (Figs. 8 to 10) to adjust the ob- 
served H/3 line width ratios by means of additional intrinsic line 
profiles and/or further line broadening mechanisms. Figure 8 
shows the computed trend of intrinsic Gaussian profiles that are 
broadened by rotation. The corresponding line profiles are pre- 
sented in Fig. 5. The intrinsic Gaussian profiles had line widths 
(FWHM) ranging from 500 to 5,000 km s~ l (from top to bottom). 
The rotational velocities go from 500 to 7,000 km It is evi- 



dent that these computed H/3 line width ratios cannot explain the 
observations because these broadened Gaussian profiles always 
exceed FWHM/cr;,„ e ratio values of 2.35. Lorentzian profiles, on 
the other hand, can easily explain low FWHM/<x tolt , ratios as seen 
in Fig. 7. 

Further line broadening tests of Lorentzian profiles con- 
volved with Gaussian profiles that have widths of 500 to 
10,000 km s^ 1 explain neither the high observed FWHM/cr/,- ne ra- 
tios nor the observed general trend of the FWHM/<x/ inf , to line 
width FWHM ratios (Fig. 9). The corresponding line profiles are 
given in Fig. 6. 

Figure 10 shows the trend of both observed and modeled 
H/3 line width ratios FWHM/cr;„, e versus FWHM for logarithmic 
profiles. We modeled logarithmic profiles for different v mi „/v max 
ratios (0.01 - 0.9, from bottom to top). In addition we varied 
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2500 kms Lorentz 



2000 4000 6000 8000 
FWHM [kms -1 ] 



10000 



12000 



x 




1000 



3000 



5000 



Fig. 11. Observed Hy line width ratios FWHM/cr/, ne versus line Fig. 14. Observed and modeled C in] /11909 line width ratios 
width FWHM and their modeling by rotational broadening of FWHM/cr ; ,„ e versus line width FWHM. 
Lorentzian profiles. Less reliable measurements are marked by 
open red diamonds. 
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FWHM [kms" 1 ] Fig. 15. Observed and modeled Si iv ,11400 line width ratios 
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the line widths from 300, 500, 700, 900, to 1 100 km s -1 (left to 
right). The corresponding line profiles are given in Fig. 2. 

An exponential profile has a similar shape to a Lorentzian 
profile, especially in their line wings (see Fig. 1). Line broad- 
ening due to rotation of these two profile types leads to similar 
trends in the FWHM/cr/,„ e vs. FWHM figures. Exponential pro- 
files have a fixed FWHM/o-, ; „ t ratio of FWHM/cr,,„ e = V2 In 2 * 
0.9803, while this ratio goes to zero for Lorentzian profiles. All 
in all, the observed trend of varying line shape as a function of 
line width - from the narrowest to the broadest line profiles - 
cannot be explained by one varying parameter alone. 

3.2.3. Rotational line broadening of Lorentzian Hy, Ha, Lya, 
C in] ,11909, He ii ,11640, and Si iv ,11400 profiles. 

All observed rms emission line width ratios FWHM/o7,„ t , versus 
line width FWHM (from Table 1), as well as their modeling, are 
shown separately in Figs. 1 1 to 16 for the different emission lines 
Hy, Ha, Lya, C m] ,11909, He n ,11640, and Si iv ,11400. The ob- 
served data are taken from the AGN sample (see Peterson et al., 
12.004b . Those measurements that were regarded as less reliable 
by the authors are marked by open red diamonds. These less 
reliable measurements follow in most cases the general trend. 
Besides the least square fits we show lower and upper limits to 
the widths of the Lorentzian profiles (in Figs. 11 to 16). These 
limits are calculated from the variation in the individual data 
points with respect to the least square fit. 

There are two Ha measurements not following the trend: the 
galaxies NGC 4593 and PG 1426+015. The galaxy NGC 4593 
shows a high line width ratio FWHM/o7, Jlf of 2.7. This might 
be explained with a high-inclination angle of the accretion 
disk (see Paper I) in this galaxy. PG 1426+015 shows a low 
FWHM/cr/„, e line width ratio with respect to the derived line 
width of 5450 km s . The same galaxy was an outlier in the cor- 
responding H/3 figure (see Paper I). 

Furthermore, there is one Si iv ,11400 measurement (Fig. 15) 
that does not follow the general trend by showing a high 
FWHM/cr;„ IC line width ratio. However, the error of this line 
measurement in the galaxy NGC 5548 is very large. A second 
measurement of this line in NGC 5548 at FWHM/o7,„ e = 1.755 
(see Table 1) corresponds with the measurements of other galax- 
ies. 

The line width ratios FWHM/cr/, ;lf , versus line width FWHM 
of all these optical and UV emission lines can be modeled in 
a simple way by rotational line broadening of Lorentzian pro- 
files, as presented before in Paper I for the H/3, Hen ,14686, 
and Civ ,11550 lines. To each emission line belongs a dedi- 
cated turbulent velocity that can be derived from the underlying 
Lorentzian profile. This specific turbulent velocity of the indi- 
vidual lines is the same in all galaxies. The intrinsic turbulent 
velocities belonging to the individual emission lines are listed in 
Table 2. 

In Paper I we considered the possible effect of an inclination 
of the accretion disk. An inclination leads to smaller line widths 
in comparison to the intrinsic line widths. Here we did not cor- 
rect for this in Table 2. 

4. Discussion 

4.1. AGN line profiles 

In early studies of emission line profiles e.g. Blumenthal & 
Matthews J1975I ) or Capriotti et al. dl980b fitted logarithmic pro- 
files to observed quasar emission lines. However, later on it was 



Table 2. Intrinsic turbulent velocities connected to the line emit- 
ting regions of the strongest emission lines. 



emission line turbulent velocity 

[kms- 1 ] 



H/3 


400 (- 


■175) (+300) 


Hy 


425 (- 


■100) (+125) 


Ha 


700 (- 


■300) (+600) 


Hen ,14686 


900 (- 


■300) (+200) 


Cm] ,11909 


1,500 (- 


■600) (+1,000) 


Si iv A 1400 


2, 100 (- 


■800) (+1,600 


Hen,U640 


2,300 (- 


■800) (+1.200) 


Civ ,11549 


2,900 (- 


■1,000) (+1,400) 


Lya+Nv ,11240 


3,800 (- 


•1,100) (+2,000) 



noticed that their logarithmic fits result in line wings that are too 
small in comparison to the observed profiles. In those early years 
only very few high-quality AGN spectra existed. 

Afterwards it became apparent that the profiles of narrow- 
line Seyfert 1 galaxies with H/3 line widths (FWHM) of less 
than 2000 - 4000 km s are well adjusted by single Lorentzian 
profiles (Sulentic et al. [2000] and references therein, Veron et 
al. 120011 Marziani et al. :20Q3). These authors confirmed previ- 
ous claims that observed broad Balmer line profiles are more 
precisely fitted by Lorentzian rather than Gaussian profiles. 
Gaussian profiles cannot reproduce the profiles observed in 
narrow-line Seyfert 1 galaxies. 

Laor d2006t found evidence of exponential line wings in 
the Ha line of the low luminosity narrow-line Seyfert 1 galaxy 
NGC 4395. However, this galaxy is unique in having a broad Ha 
line width (FWHM) of < 520 km s _1 only. Laor ( I20U61 ) explains 
this specific line profile by electron scattering in the broad-line 
region. 

Observed H/3 profiles in AGN with line widths (FWHM) < 
4,000 km s~ l never can be explained by single Gaussian profiles 
or by a combination of Gaussian profiles. Gaussian profiles have 
a constant FWrlM/crj,,, ratio of 2.35. In contrast, all observed 
H/3 profiles showing line widths (FWHM) < 4,000 kms~ l ex- 
hibit lower FWHM/<x to , ? ratios (see Fig. 8). A combination of 
Lorentzian profiles with Gaussian profiles does not match the 
observed trend in the FWHM/cr/,, I( , versus line width FWHM 
plots (see Fig. 9). Logarithmic profiles meet the observed trend 
in the FWHM/<x/,„ P versus FWHM plots for those line profiles 
with line widths FWHM < 6,000 kms~ l (Fig. 10). However, 
broader profiles cannot be modeled in a simple way by varying 
only one parameter. 

The profiles are not the same in all the quasars. Sulentic et 
al. d2000b call those quasars having H/3 line widths (FWHM) > 
4,000 kms~ l Population B quasars in contrast to Population A 
quasars having H/3 line wi dths (F WHM) < 4,000 kms~ l . In a 
similar spirit, Collin et al. (120061 call those galaxies emitting 
line profiles narrower than Gaussian profiles ( FWHM/cr/, n( , < 
2.35) Population 1 AGN and those galaxies showing line profiles 
broader than Gaussian profiles Population 2 AGN (their Fig. 3). 
Sulentic et al. (2000) modeled the broad H/3 line profiles (i.e. the 
Population B quasars) with two Gaussian profiles, a broad one 
and a narrower one. In no case could single Gaussian profiles 
match their observed broad-line wings. 

We demonstrated in Paper I that the shape (i.e. the 
FWHM/cr /ine ratio) of the H/J, Hen ,14686, and Civ ,11550 pro- 
files varies systematically with their line width. There is no clear 
evidence for two separate quasar populations. Here we show 
that this applies for Hy, Ha, Lya, He n ,11640, Cm] ,11909, and 
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S iv/11400 lines in our sample as well. The observed profiles of 
all these emission lines can simply be characterized by rotational 
line broadening of Lorentzian profiles. To every emission line 
belongs one exclusive turbulent velocity of the line emitting re- 
gion. The turbulent velocities go from 400 km s _1 for the low- 
ionization lines up to 3,800 km s~ l for the high-ionization lines 
(see Table 2). The rotational velocities go from 500 km s up to 
6,500 fern*- 1 . 

Baskin & Laor (2005) claim that narrow CIV lines are rare 
(~2 per cent occurrence rate) compared with narrow H/3 < 
2,000 km s~ l (~20 per cent) based on more than 80 spectra from 
the Boroson-Green sample. This can be understood easily in 
terms of the high turbulent velocities belonging to the C iv /il550 
line emission region (2,900 km s~ l ) in comparison to the H/3 line 
region (400 km s~ l ). 

Line broadening due to turbulence and rotation are the 
main constituents for the observed broad-line profiles in AGN. 
However, inclined accretion disk geometries of the line emit- 
ting regions lead to smaller line widths owing to projection 
effects, while their FWHM/cr; !nc ratio remains constant (see 
Paper I). Additional asymmetries in the line profiles might be 
caused by geometrical/optical obscuration effects, additional 
outflow/inflow components, anisotropic emission, superposi- 
tion of line emission from different emitting region, etc. It 
was mentioned in Paper I that line asymmetries lead to lower 
FWHM/crune ratios as well. 

This general topic should be studied in more detail in the 
future. Furthermore, obscuration affects individual line profiles 
in a different way because the lines originate in different regions 
(see below). 

4.2. Geometry and structure of the line emitting region 

In Table 2 one can identify the clear trend for higher ionized 
lines to originate in those regions where higher turbulent veloc- 
ities are predominant. As higher ionized lines exhibit broader 
emission lines in general and usually originate closer to the cen- 
ter (e.g. Peterson & Wandel [l999l Kollatschny et al., l200TT >. This 
trend is consistent with a general increase in the turbulent veloc- 
ity towards central regions. The broad emission lines in AGN 
originate at distances of less than one light day to more than 100 
light days from the central ionizing source (e.g. Desroches et al. 
[20061 Kaspi et al. [20071 Bentz et al. |2009l . 

Based on the earlier theoretical studies of Pringle (119811 Eq. 
3.16), we made the following claim about the geometry of AGN 
accretion disks in Paper I. The ratio of the accretion disk height 
H with respect to their radius R is proportional to the turbulence 
velocity v tU rb m the accretion disk with respect to the rotational 
velocity v rot 

H/R = (l/a)(v turb /Vrot). 0) 

The unknown viscosity parameter a is assumed to be con- 
stant. We did not consider magnetic fields in this picture and 
made the obvious assumption that the turbulence velocity v tur b is 
less than the sound of speed c$ in the disk: 

V,urb < C S (4) 

v lurb = a' c s (5) 

where a' is close to the usual a parameter. 

The rotational velocities - belonging to the individual emis- 
sion lines - vary by a factor of more than ten, while the turbu- 
lent velocities that are connected with the individual emission 
lines remain constant. This can be seen in Fig. 7, as well as in 



Fig. 17. Schematic accretion disk models of AGN showing 
that slow-rotating AGN have an accretion disk thst is ten times 
thicker than fast-rotating AGN. 

Figs. 1 1 to 16. In this accretion disk model and accepting our re- 
sults breaking the velocity degeneracies one would conclude that 
slow-rotating AGN host a thick accretion disk that is ten times 
thicker than fast-rotating AGN. A schematic picture of thick and 
thin accretion disks is shown in Fig. 17. The radius of the cen- 
tral black hole (r = 5.9 x I0 n cm) in Fig. 17 corresponds to a 
Schwarzschild mass of M = 2 x 10 7 M o . 

Other details regarding the physical conditions in the line 
emitting region seem to be even more complex: It has been 
noticed by means of reverberation measurements that the 
He H/14686 and Hen A 1640 lines in AGN spectra originate at 
different distances from the center (e.g. Peterson et al., I2004t . 
although they hold the same ionization degree. Diverse mod- 
els do not reproduce this observational fact (e.g. Bottorff l2002l >. 
Now we independently confirm the former finding that these 
lines originate in different physical regions based on the dis- 
tinct turbulent velocities we deduced for the two helium II 
line emission regions. This is important considering that the 
He ii /14686/He n /1 1640 ratio has sometimes been used as a red- 
dening indicator for the broad-line region (e.g. Snijders et al. 
[T9861 Ferguson et al. [19951 

The eigenvector studies of Boroson & Green (1992) demon- 
strated a strong correlation between the emission line width, 
optical Fell emission, and soft X-ray photon index in AGN 
(Eigenvector 1). The line width is connected with the geometry 
or rather, to be more specific, with the thickness of the accretion 
disk (see above). Therefore, the observed strength of the spectral 
Fell emission and the X-ray photon index might be affected by 
geometry effects as well. 

4.3. Correction factors for calculating central black hole 
masses 

The central black hole mass Mbh in AGN can be derived from 
the broad emission line widths under the assumption that the gas 
dynamics are dominated by the central massive object: 

M BH =fcT cenl Av z G- 1 . (6) 

The characteristic distances of the line-emitting regions 
c T ce „, from the central ionizing source can be derived by means 
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of spectroscopic variability campaigns. The distance can be cal- 
culated from the cross-correlation function of emission line in- 
tensity variations with respect to the ionizing continuum inten- 
sity variations (e.g. Koratkar & Gaskell 119911 Kollatschny & 
Dietrich 1997). The characteristic velocity Av of the emission 
line region can be estimated from the FWHM of the rms profile 
or from the line dispersion cr/,„ e (e.g. Peterson et al. [2004 ). The 
scale factor f depends among others on the geometry and struc- 
ture of the line emitting region, as well as on their inclination. 
Scale factors for rms spectra have numerical values from 0.5 up 
to 6.2 (e.g. Collin et al. 12006k depending on the assumptions 
about geometry and structure of the line emitting region. Only 
very few cases are known where the black hole mass could be 
estimated independently by, e.g., gravitational redshift measure- 
ments (Kollatschny 2003b) to get additional information on the 
scale factor f or rather the inclination angle. 

Variability campaigns are very expensive with respect to 
their observing time. Kaspi et al. (120051 1 verify a relationship 
between AGN luminosity and Balmer-line-averaged BLR size 
based on variability campaigns. This scaling relation can be used 
to obtain black hole masses based on single epoch spectra alone. 
The velocity of the line emitting region can be estimated on the 
basis of the width of the Balmer emission lines. This method was 
extended to other emission lines, especially to the Civ A 1549 
line of distant high-redshift AGN. 

Central black hole masses were generally found to be bigger 
for AGN showing broad H/3 profiles (Pop. B), as well as broad 
Crv/(1550 profiles than for those AGN that exhibit narrow H/3 
profiles (Po p. A) (Marziani et al. [20031 Peterson et al., 120041 
Vestergaard l2004l Zamfir et al. l20T0l 

The broader emission lines usually originate closer towards 
the center as has been noted before (e.g. Peterson et al. 120041 
Kollatschny 12003 al ). The turbulence grows towards the inner 
zones as well. Therefore the black hole mass estimates, based 
on the line widths, are more heavily biased in broad emission 
line objects caused by the additional turbulence broadening than 
in the narrower emission line objects. 

We consider here nearby AGN where we have good black 
hole mass estimations from reverberation mapping. Shen et al. 
( 2008 ) have investigated biases in virial black hole masses based 
on SDSS spectra. Wang et al. ( 1201 11 1 discuss the influence of 
radiation-driven outflows on CIV emission line profiles in high- 
redshift and high-luminosity AGNs. In the case of such outflows, 
the CIV emission lines should become broader because of the 
additional line broadening component of the outflow. If these 
CIV profiles of high-redshift and high-luminosity AGNs show 
the same widths as those of nearby AGNs, then the contribution 
of the rotation on the line width must be even lower. In that case 
the black hole mass estimation that is only based on the rotation 
velocity must be scaled down even more. 

Vestergaard d2004t derived black hole masses of nearby 
quasars based on their H/3 line widths, as well as black hole 
masses of high-redshift quasars based on their Civ A 1549 line 
widths. On average, she found bigger black hole masses - by a 
factor five to ten - for the younger high-redshift quasars based 
on their C iv /11549 line widths in comparison to the old nearby 
quasars based on their H/3 line widths. In an additional paper, 
Vestergaard & Peterson ( 2006) published different scaling rela- 
tions for nearby and distant AGN black hole masses based on 
their H/3, as well as Civ Al 549 line widths. However, Netzer at 
el. (I2007t note that using the Civ A 1549 line width for estimat- 
ing black hole masses gives considerably different results and a 
larger scatter than using the H/3 line. We demonstrated in Paper I 
that black hole mass estimates based on the line width of the 
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line emitting region from the observed FWHM, shown for ob- 
served profiles having intrinsic turbulence velocities ranging 
from 400 km s" 1 up to 3,800 km s" 1 . 



Civ/11549 line are overestimated by a factor of five to ten in 
comparison to those black hole masses based on the width of the 
H/3 line. 

On the basis of our observed and modeled line width ratios 
in Figs. 7, 11 - 16, we present in Fig. 18 the relation between 
the observed FWHM of the emission line profiles and the re- 
lated rotational velocity v rot . This relation is shown for under- 
lying turbulence profiles with v„„;, ranging from 400 kms~ l to 
3,800 km s" 1 . In Fig. 19 we present the correction factor towards 
the rotational velocity v rot of the broad-line clouds to calculate 
the central black hole masses of AGNs based on the observed 
line widths. 

Computed black hole masses go with the square of the ro- 
tational velocity (see Eq. |6j. Therefore, the corrected intrinsic 
black hole masses should be lower by a factor of two to more 
than ten with respect to the calculated black hole masses that 
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have not been corrected for the turbulent velocity contribution. 
This correction factor is different for the different emission lines 
(depending on the underlying turbulent velocity of the individ- 
ual lines) and also depends on the additional rotational velocity 
broadening. There is not an exclusive correction factor for black 
hole mass estimates based on the C iv /U549 or H/3 line widths as 
proposed by Vestergaard & Peterson (120061 in their mass scaling 
relations. Furthermore, black hole masses of distant AGN based 
on the CiV/11549 line widths have been overestimated in com- 
parison to nearby AGN where the mass estimate is based on the 
H/3 line. It is important to consider this effect for our general un- 
derstanding of the evolutionary history of black hole masses in 
AGN. 

5. Conclusion 

We investigated the profile shapes of the UV/optical broad emis- 
sion lines in AGN in detail. The two basic components caus- 
ing the line profile are Lorentzian profiles and rotational broad- 
ening. An intrinsic turbulent velocity belongs to each specific 
AGN emission line, which causes different line widths (FWHM) 
of the particular Lorentzian profiles. The turbulent velocities go 
from 400 kms' 1 for H/3 up to 3,800 km s" 1 for Lycr+Nv ,11240. 
The rotation velocities causing the line profile broadening go 
from 500 km s~ l up to 6,500 km s . There are not two separated 
classes of broad-line AGN (narrow line and broad-line objects), 
but instead a continuous transition from narrow to broad-line ob- 
jects. The accretion disk thickness in AGN can be derived from 
the ratio of the turbulent velocity v tU rb with respect to the rota- 
tional velocity v rot . Slow-rotating AGN have an accretion disk 
that is ten times thicker than fast-rotating AGN. We find clear 
evidence that the Hen/14686 and Hen A 1640 lines originate in 
different physical regions in AGNs, although they hold the same 
ionization degree. This finding is based on the distinct turbulent 
velocities we deduce for the two helium II line emission regions. 
Black hole mass corrections resulting from turbulent velocity 
considerations are nontrivial and are larger for C iv /11549-based 
measurements than H/3, and both are a function of the observed 
FWHM. In the literature usually the widths of the broad emis- 
sion lines are used to compute the central black hole masses in 
AGN. However, one has to consider the contribution by the tur- 
bulence to get the corrected intrinsic rotational velocities. This 
investigation presented the individual correction factors towards 
the different emission lines for getting the intrinsic rotational 
velocities v rot of the line emitting region from the observed 
FWHMs. The corrected black hole masses are lower than the 
uncorrected black hole masses by a factor of two to more than 
ten, depending on the emission lines, as well as on the rotational 
velocities. Especially those masses of the distant AGN have been 
grossly overestimated where the masses have been estimated on 
the basis of the C iv /11549 line widths. 
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